Phase diagram of dense neutral three-flavor quark matter 
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We study the phase diagram of dense, locally neutral three-flavor quark matter as a function of 
the strange quark mass, the quark chemical potential, and the temperature, employing a general 
nine-parameter ansatz for the gap matrix. At zero temperature and small values of the strange 
quark mass, the ground state of matter corresponds to the color-flavor-locked (CFL) phase. At 
some critical value of the strange quark mass, this is replaced by the recently proposed gapless CFL 
(gCFL) phase. We also find several other phases, for instance, a metallic CFL (mCFL) phase, a 
so-called uSC phase where all colors of up quarks are paired, as well as the standard two-flavor 
color-superconducting (2SC) phase and the gapless 2SC (g2SC) phase. 



I. INTRODUCTION 

At sufficiently high densities and sufficiently low tem- 
peratures quark matter is a color superconductor [J] . This 
conclusion follows naturally from arguments similar to 
those employed in the case of ordinary low-temperature 
superconductivity in metals and alloys Of course, 
the case of quark matter is more complicated because 
quarks, unlike electrons, come in various flavors (e.g., 
up, down, and strange) and carry non-Abelian color 
charges. In the last six years, this phenomenon was stu- 
died in detail by various authors [all H Irj U.M.WL Ho| . 
Many different phases were discovered, and recent studies 
0,0,0,0,0,0] suggest that even more new phases 
may exist. (For reviews on color superconductivity see, 
for example, Ref. 0].) 

In nature, the most likely place where color super- 
conductivity may occur is the interior of compact stars. 
Therefore, it is of great importance to study the phases 
of dense matter under the conditions that are typical 
for the interior of stars. For example, one should appre- 
ciate that matter in the bulk of a star is neutral and f3- 
cquilibratcd. By making use of rather general arguments, 
it was suggested in Ref. [2lJ that such conditions favor 
the color-flavor-locked (CFL) phase and disfavor the two- 
flavor superconducting (2SC) phase. In trying to refine 
the validity of this conclusion, it was recently realized 
that, depending on the value of the constituent (medium 
modified) strange quark mass, the ground state of neutral 
and /3-equilibrated dense quark matter may be different 
from the CFL phase 0. In particu- 
lar, the _gapless two-flavor color-superconducting (g2SC) 
phase 0| is likely to be the ground state in the case of 
a large strange quark mass. On the other hand, in the 
case of a moderately large strange quark mass, the reg- 
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ular and gapless color-flavor-locked (gCFL) phases 0] 
are favored. At nonzero temperature, T 7^ 0, some other 
phases were proposed as well 0. 

In this paper, we study the phase diagram of dense 
neutral three-flavor quark matter as a function of the 
strange quark mass at zero and finite temperature. In or- 
der to allow for the most general ground state, we employ 
a nine-parameter ansatz for the gap function. The effects 
of the strange quark mass are incorporated in our model 
by a shift of the chemical potential of strange quarks, 
Ms ~ *■ Ms — m s/(2M) where i = r, g,b is the color index, 
m s is the strange quark mass, and y, is the quark chemical 
potential. This shift reflects the reduction of the Fermi 
momenta of strange quarks due to their mass. Such an ap- 
proach is certainly reliable at small values of the strange 
quark mass. We assume that it is also qualitatively cor- 
rect at large values of the strange quark mass. 

We note that the analysis of this paper is restricted to 
locally neutral phases only. This automatically excludes, 
for example, mixed [Isj and crystalline |l9j | phases. Also, 
in the mean field approximation utilized here, we cannot 
get any phases with meson condensates [2pj . 

This paper is organized as follows. In Sec. [H] we 
present the details of our approach which is based on 
the Cornwall- Jackiw-Tomboulis (CJT) formalism [22|| . 
There, we also derive the gap equations and the neu- 
trality conditions, and obtain an expression for the pres- 
sure. In the next two sections, the gap equations and 
the neutrality conditions are studied by using numerical 
methods. The results at zero temperature are presented 
and discussed in Sec. IIIII There, the appearance of the 
CFL and gCFL phases is established in the case of small 
and moderately large values of the strange quark mass, 
respectively. The results at nonzero temperature are dis- 
cussed in Sec. II VI Also, the phase diagram of quark mat- 
ter in the T-m s / [i plane, as well as in the T-\i plane 
is presented. Finally, Sec. IVl concludes this paper with a 
summary of the results. 

Our units are % = c = ks = 1. The metric tensor 
is g^v — diag (1, — 1, — 1, — 1). Four- vectors are denoted 
by capital Latin letters, e.g., — (fc°,k) where k is a 
three- vector with absolute value k = |k| and direction 
k = k/fc. We use the imaginary-time formalism, i.e., the 
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space-time integration is denned as J x = J^ T dr J v d 3 x, 
where r is the Euclidean time coordinate and V the three- 
volume of the system. Energy-momentum sums are de- 
fined as follows: T/VJ2k = T E„J d 3 k/(27r) 3 where 
the sum runs over the fermionic Matsubara frequencies 
w n = (2n + l)irT = ik(). 



II. MODEL AND FORMALISM 

In this paper, we follow closely the approach of Ref. [2^ 
which is based on the C JT formalism [22] • We generalize 
the model of Ref. in the way that strange quarks may 
now also participate in forming Cooper pairs. The quark 
spinor field has the following color-flavor structure: 

v. 
r u 

v> = r d . (i) 

The Dirac conjugate spinor is defined as t/> = V> 7o- In the 
treatment of superconducting systems it is advantageous 
to double the fermionic degrees of freedom by introducing 
the so-called Nambu-Gorkov spinors 



* = 



(2) 



where ipc = Cip T is the charge-conjugate spinor, and C 
is the charge-conjugation matrix. 

In this paper, as in Ref. |23| . we shall approximate 
the gluon-exchange interaction between the quarks by a 
point-like current-current interaction. Physically, this is 
equivalent to a model with heavy, non-dynamical gluons. 
In many ways, such an approximation is analogous to the 
nonrenormalizable Fermi theory of weak interactions that 
was used before the SU(2)xU[l) gauge theory of electro- 
weak forces was developed 24]. Also, in the mean-field 
approximation that we use below, our treatment will be 
similar to that in Refs. @ S [H El III 111 E3 
which are based on Nambu-Jona-Lasinio (NJL) [28| type 
models. 

In the approximation with a point-like interaction, the 
CJT effective action simplifies. It has the following gene- 
ral form: 



f 



1 



r [S] = -TrlnS- 1 + -Tr (S^S - l) +T 2 [S] 



(3) 



Here the traces run over space-time, Nambu-Gorkov, 
color, flavor, and Dirac indices. The factor 1/2 in front of 
the fermionic one-loop terms compensates the doubling 
of the degrees of freedom in the Nambu-Gorkov basis. 



In the above effective action, S denotes the full quark 
propagator. 

In order to understand the structure of the full quark 
propagator S, let us first discuss the structure of the free 
quark propagator So in the Nambu-Gorkov basis. The 
inverse, S^ 1 , is given by the following matrix: 







where 



[G - 



(4) 



(5) 



are the inverse Dirac propagators for massless quarks 
and charge-conjugate quarks, respectively. At sufficiently 
large quark chemical potential, there is no need to take 
into account the small masses of the light up- and down- 
quarks. It is easy to understand that the dynamical effect 
of such masses around the quark Fermi surfaces is negli- 
gible. Of course, the situation with the strange quark is 
different because its mass is not very small as compared 
to /i. It appears, however, that the most important ef- 
fect of a nonzero strange quark mass is a shift of the 
strange quark chemical potential due to the reduction 
of the Fermi momentum, /i* — * fi\ — m 2 s /{2[i). (Strictly 
speaking, it is n\ rather than \i that should appear in 
the denominator. Quantitatively, however, this does not 
make a big difference.) For simplicity, this is the only 
mass effect that we include in our analysis below. Note 
that nonzero quark masses were properly accounted for 
in Refs. HUHE!. 

The quark chemical potential matrix jl in color-flavor 
space is defined as 



p, = diag 



(6) 



In realistic systems, such as bulk matter inside a star, 
some of the ify should be related. This is because, in 
chemical equilibrium, one can introduce only as many 
independent chemical potentials as there are different 
conserved charges in the system. In quark matter, the 
chemical potentials im\ can be defined in terms of the 
quark chemical potential fi (fi = /is/3, where \ib is the 
baryon chemical potential) , the chemical potential fiQ for 
the electrical charge, and the two chemical potentials, /Z3 
and fig, for color charge, 



H) = fx + n Q Q f + fi 3 T£ + fi 8 Ts 



(7) 



Thus, there are only four out of nine chemical potentials 
in Eq. © that are independent in the case of dense three- 
flavor quark matter in chemical equilibrium. In (locally) 
neutral matter, as we shall see below, only one of them 
will remain independent. 

In the framework of the CJT formalism, the thermo- 
dynamic potential of quark matter is proportional to the 
CJT effective action at its stationary point, determined 
by the solution of the following equation: 



ST 

ss 



0. 



(8) 
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This is nothing but the Dyson- Schwinger equation for the 
quark propagator, 



where 



— 5n 



[S] 



£ = 2 



6S 



(9) 



(10) 



is the quark self-energy. The functional r 2 is the sum 
of all two-particle irreducible (2PI) diagrams. Unfortu- 
nately, it is impossible to evaluate all 2PI diagrams ex- 
actly. Nevertheless, the advantage of the CJT effective 
action @ is that a truncation of the sum r 2 at a finite 
number of terms still provides a well-defined many-body 
approximation. In this study, we only include the sunset- 
type diagram shown in Fig. 1 of Ref. [23|. which becomes 
a double-bubble diagram in the case of a local, instanta- 
neous interaction. It is easy to check that this leads to 
the following expression for the self-energy: 



E(#) = -s 2 ^r^(Q)r^ 



v 



(ii) 



The regular (diagonal) part of the self-energy is neglected 
here. While it plays an important role in the dynamics 
of chiral symmetry breaking, it is essentially irrelevant 
for color superconductivity in dense quark matter. (The 
effect of the diagonal part of the self-energy was studied, 
for example, in Ref. The off-diagonal elements <f>~ 
and $ + in Eq. i|14|) are matrices in Dirac, color and fla- 
vor space. They are related by the requirement that the 
action is real, $~ = 7 ($+) t 7 . 

With Eq. (jnj one may invert Eq. @ to obtain the full 
quark propagator, 



S 



G+ 5- 
S+ G- 



(15) 



where the diagonal and off-diagonal elements are 



and 



G ± = {[G ± ]- 1 -^ G ^ ± }- 1 



(16) 



Local, instantaneous gluon exchange is parametrized by 
a propagator of the form 



7-jab _ cab 9/J.u 

— A 2 ' 



(12) 



In Eq. I|ll|) . we introduced the Nambu-Gorkov vertex, 

(13) 



l»T a 
-fT? 



Equation Ijlljl is a self-consistency equation for £, com- 
monly called the "gap equation" . 

In color-superconducting systems, the quark self- 
energy has the form 





$+ 



(14) 



(17) 



respectively. 



In the following, we use the notation T* for the CJT 
effective action evaluated at the stationary point. In the 
approximation used here, it takes the following form: 



T* = -TrlnS 
2 



1 ■ -Tr CSS) 
4 



(18) 



In this paper, we utilize the following nine-parameter 
ansatz for the gap matrix: 
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[Af s ] ± / 



(19) 



The explicit Dirac structure of the matrix elements is 
1 " ( h ) ■■ 

c,e— ± 



(20a) 



[A&,] (K) = Y €?f> {K)VZ e M ■ (20b) 



c.e— ± 



In this representation, 4ffip are real-valued gap para- 
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meters, while 



7> c e (k) = i(l + c 75 )(l + e7o7-k) 



(21) 



are the energy- chirality projectors. 

In order to calculate the first part of the effective ac- 
tion in Eq. l|18fl . we transform the inverse quark propa- 
gator © into a block-diagonal form which is schemati- 
cally shown in Fig. ^ This is easily achieved by chang- 
ing the order of rows and columns in color, flavor, 
and Nambu-Gorkov space. The energy-chirality projec- 
tors simplify the calculation of the Dirac traces. Af- 
ter performing these traces, we use the matrix relation 
Tr In A = J2k m (^ e * A) *° calculate the one-loop contri- 
bution of the quarks to T*. 



6x6 



2x2 



2x2 



2x2 



2x2 



FIG. 1: The block-diagonal structure of the inverse quark 
propagator in color-flavor and Nambu-Gorkov space. 



We use the Gauss elimination procedure to reduce the 
6 x 6-block to a 3 x 3-block and all the 2 x 2-blocks to 
1 x 1-blocks in color-flavor space. The determinant of 
the 3 x 3-block is computed analytically by using, for 
example, Mathematica or Maple. 

The most complicated expression arises from the de- 
terminant of the 3 x 3-block. Schematically, it has the 
following form: 



det (3 x 3) e = k% + bk% + cfcg + d. 



(22) 



Here, the coefficients b, c, and d are rather compli- 
cated functions of the quark momentum k, three che- 



mical potentials (/^, n 9 d , and /i h s ) and six gap parame- 
ters f u % and <t>f s ). This determi- 

nant can always be factorized as follows: det (3 x 3) e = 

[*o - (efTPo - (£2) 2 Po - (el) 2 ]- A « k cl ear, the func- 
tions if, l\ and e| determine the dispersion relations of 
the three quasiparticles described by the 3 x 3-block of 
the quark propagator. In order to get their explicit ex- 
pressions, one solves the cubic equation det (3 x 3) e = 
for £ = fcg. By making use of Cardano's formula, the so- 
lutions can be presented in an analytical form that we 
use later in the numerical calculations. Because of their 
very complicated nature we refrain from presenting them 
explicitly. 

The dispersion relations of the quasiparticles that cor- 
respond to the 1 x 1-blocks are easy to derive. They 
are given in terms of the same function, e£(/x, (ft) — 
[(k — efi) 2 + l^l 2 ] 1 / 2 , that also appears in the 2SC phase. 

Once the quasiparticle dispersion relations are estab- 
lished, the calculation of the first term in the CJT action, 
i.e., J2k ln (detS- 1 ) = £\ £ K In [(w n + ifij) 2 + e 2 ] , re- 
duces to a sum of several standard contributions of the 
following type |3(| : 



K 



j"2 



V 



p 2 dp 
~2~^~ 



T 



+ ln (l + e -fo-«>/ r ) + ln (l + e -< e '+'*>/ T ) 



(23) 



In order to simplify the second term in the effective po- 
tential in Eq. I)18[l. we use the gap equation (| 1 1 1> . Then, 
we arrive at the following simple result: 



-iTr(E5) = ^^Tr 
4 1 ' 4g 2 T 



(24) 

where the new notation $ was introduced to represent 
only the color- flavor part of the gap matrix <i> + , given 
in Eq. (|19fl . In this connection, it should be noted that 

the gap parameters 0c//' are independent of the chirality 
and energy projections. This is a consequence of the gap 
equation i|ll|) . 

Finally, by combining all contributions to Eq. <|18[) . we 
derive the result for the pressure p = pT*, 
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T sr 

i 



7T- 



dk k 2 I In 



1 + exp 



Ei - fM 
T 



In 



1 + exp 



Ei + m 

T 



3 A* 
4? 



E 



i=le=± Jo 1 



1 + exp 



= 1 e=± 



^EE ( dfefe a ^e£Oi i ,^)-* + riii 

+Tln 



l + e *p +^ 



1 + exp - 



T 
T 



(25) 



The first term in this expression is the contribution of lep- 
tons (i.e., I — e~,(jT). In principle, the contribution of 
neutrinos should be added as well. In this paper, however, 
their contribution is neglected. This is a good approxima- 
tion for compact stars after deleptonization. The disper- 
sion relations of the leptons are given by Ei = ^Jk 2 + mf 
with m e - w 0.511 MeV and m„- « 105.66 MeV. Note 
that the vacuum contribution was subtracted in Eq. (|25H , 
in order to improve the convergence of the momentum in- 
tegrals. In various applications, a bag constant could be 
added if necessary. In Eq. (25)1 , <f> t are the three indepen- 
dent gap parameters that appear in the 2 x 2-blocks in 
Fig. ^ while ipi and <7j appear in the 6 x 6-block. They 
are defined by 

0i = <t% b d> = 0ds. and °"i = 0««> (26a) 

<h = <f&, = and a 2 = <j> 9 d % (26b) 

03 = 4>2,, ¥>3 = f u % and °3 = 4>f s - (26c) 

In accordance with the fact that the gap parameters 
are independent of the energy and chirality, the corre- 
sponding indices were dropped. Similarly, jii and Sfii 
(i = 1,2,3) are the average values and the differences 
of various pairs of chemical potentials that come from 
the three different 2 x 2-blocks, 

fii = -(m! + and 6m = -(fi 3 s - fi b d ), (27a) 

M2 = \(fi r a + ft), and 5^2 = \(n r a - ft), (27b) 

= + and Sfj,3 = -(fid - (27c) 

In order to render the integrals in the expression for the 
pressure finite, we introduced a three-momentum cut- 
off k. In QCD with dynamical gluons, of course, such 
a cutoff would not be necessary. Here, however, we use 
a model with a local current-current interaction which is 
nonrenormalizable . 

In the case of the ansatz (I19|) for the gap function, there 
are nine coupled gap equations (corresponding to the to- 
tal number of gap parameters) that follow from Eq. (jllfl . 
The explicit form of these equations is very complicated 



and not very informative by itself. Because of this, they 
will not be displayed here. In the next two sections, we 
solve these gap equations by numerical methods. 

As was mentioned in the Introduction, matter in the 
bulk of a compact star should satisfy the conditions of 
charge neutrality and /3-equilibrium. The latter is auto- 
matically satisfied after choosing the chemical potentials 
as in Eq. Q. The neutrality conditions read 

n Q = ^- = 0, (28a) 

OfiQ 

n 3 = p- = , (28b) 

Ofl 3 

n 8 =-^=0. (28c) 

Ofl8 

After solving these equations, one is left with only one 
chemical potential (out of four) that remains indepen- 
dent. It is most convenient to keep the quark chemical 
potential (i as a free parameter, and determine fiQ, ^3 
and /Z8 from Eqs. (j28(l . This is done numerically. 

III. RESULTS AT ZERO TEMPERATURE 

In this section we focus on three-flavor quark matter at 
zero temperature. It is clear that, for small and moderate 
values of the strange quark mass, the ground state of neu- 
tral quark matter should correspond to either the regular 
(gapped) CFL phase 3 or the gapless CFL phase [lfij . 
At very large strange quark mass and/or relatively weak 
coupling, the ground state can also be either a regular 
(gapped) or gapless 2SC color superconductor [llj |. 

Before proceeding to the results, let us specify the pa- 
rameters of the model. The strength of the diquark cou- 
pling and the value of the cutoff in the momentum inte- 
grals are fixed as follows: 

g 2 /A 2 = 45.1467 GeV~ 2 , (29a) 
K = 0.6533 GeV . (29b) 

In order to see how the phase structure of neutral three- 
flavor quark matter changes with the mass of the strange 
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quark m s , we solve a coupled set of twelve equations, 
i.e., nine gap equations and three neutrality conditions, 
given in Eq. I)28[l. for various values of m s , keeping the 
quark chemical potential fixed. In this calculation we 
take fi — 500 MeV. The results for the absolute va- 
lues of the gap parameters and the chemical potentials 
MQ) M3> and are shown in Figs. [21 and 03 respec- 
tively. Note that, strictly speaking, the gap parameters 
do not coincide with the actual values of the gaps in 
the quasiparticle spectra. In the case of the CFL phase, 
for example, there is a degenerate octet of quasiparticles 
with a gap O ctct = \4>i\ and a singlet state with a gap 
^singlet = 3ipi -\fa I . In the CFL phase, fa = fa — fa < 0, 
<Pi = ¥>2 = ¥>3 > 0, and ci = o"2 = ^3 > 0. Also, in the 
CFL phase, the following relation between the gap para- 
meters is satisfied: Oi = ifi — \fa\ = 20( 6)6 \, i = 1,2,3, 
where faj$ g) is the sextet-sextet gap in the notation of 
Ref. 0. 

The results in Figs. and extend the results of 
Ref. |15| by considering a more general ansatz for the gap 
matrix that takes into account, in particular, the pairing 
in the symmetric sextet-sextet channel. The effect of in- 
cluding pairing in the symmetric channel is a splitting 
between the pairs of gaps (\fa\, <pi), (\fa\, <p2) and (\fa\, 
<pa) that is also reflected in the change of the quasiparticle 
spectra. Of course, in agreement with the general argu- 
ments of Refs. 0, Q, the symmetric sextet-sextet gaps 
are rather small, see Fig. [U (c). This, in turn, explains 
the fact why the splittings between the above mentioned 
pairs of gap parameters are not very large [compare the 
results in Figs. El (a) and (b)]. 

In this study, we confirm that the phase transition from 
the CFL phase to the gCFL phase happens at a critical 
value of the parameter rn^j \i that is in good agreement 
with the simple estimate of Ref. , 

o 

TTi 

~ 2A - 190 MeV. (30) 

A* 

The qualitative results for the chemical potentials [iq, 113, 
and /is in Fig.|3are in agreement with the corresponding 
results obtained in Ref. 4,5] as well. In the CFL phase, 
in our notation the color chemical potential /is fulfils the 
identity fi$ — —m 2 s /{^/?,^). In addition, while the CFL 
phase requires no electrons to remain neutral, the pairing 
in the gCFL phase is distorted and a finite density of elec- 
trons appears. This is seen directly from the dependence 
of the electrical chemical potential \xq in Fig. [21 which 
becomes nonzero only in the gCFL phase. This observa- 
tion led the authors of Ref. [ljj to the conclusion that 
the phase transition between the CFL and gCFL phase 
is an insulator-metal phase transition, and that the value 
of the electron density is a convenient order parameter 
in the description of such a transition. In fact, one could 
also choose one of the differences between number densi- 
ties of mutually paired quarks as an alternative choice for 
the order parameter In either case, there does not 
seem to exist any continuous symmetry that is associated 
with such an order parameter. To complete the discussion 



of the chemical potentials, we add that the other color 
chemical potential, ^3, is zero only in the CFL phase at 
T = 0. 

The effects of a nonzero strange quark mass on the 
phase structure of neutral strange quark matter could be 
viewed from a different standpoint that, in application 
to stars, may look more natural. This is the case where 
the dependence on the quark chemical potential is stu- 
died at a fixed value of m s . The corresponding numerical 
results are shown in Fig.0] (Note once again that fa with 
i = 1, 2, 3 have negative values, and we always plot their 
absolute values.) In this particular calculation we choose 
m s = 300 MeV. 

At large values of the quark chemical potential \ji > 
to;?/ (2 A) ~ 475 MeV which is similar to the small strange 
quark mass limit considered before], the ground state of 
quark matter is the CFL phase. At smaller values of the 
chemical potential, the ground state of dense matter is 
the gCFL phase. In this case, there are nine gap parame- 
ters all of which are different from each other. One could 
also check that the density of quarks that pair are not 
equal in the gCFL phase. This can be seen from Fig. [S] 
where all nine quark number densities are plotted for 
the same value of the strange quark mass, m s = 300 
MeV. Only n r d = n g u and n r s ss n b u , all other quark num- 
ber densities are different from each other. This agrees 
with the general criterion of the appearance of gapless 
phases at T = that was proposed in Ref. ^l| m t ne 
case of two-flavor quark matter. In the ordinary CFL 
phase, in contrast, one finds that n r u = n d = n b , and 
n r d = n r s = < = n B s = n b u = n b d . 

In order to see that the gCFL phase indeed describes 
a gapless superconductor, it is necessary to show that 
the dispersion relations of quasiparticles contain gapless 
excitations. In Fig. [SJ the dispersion relations of all nine 
quasiparticles are plotted. The dispersion relations for 
the corresponding anti-particles are not shown. The dis- 
persion relations that result from the 6x6 color-flavor 
block of the quark propagator are labelled by ej in Fig. 
(see the discussion and Fig. ^ in Sec.|nJ. The remaining 
dispersion relations correspond to quasiparticles that are 
described by the 2 x 2-blocks of the quark propagator, 
given by fco = e^(p,i,fa) ± Qualitatively, these are 
the same as the dispersion relations that appear in the 
2SC/g2SC phase in Ref. [HI- 

This is not a coincidence 
because the corresponding 2x2 color-flavor blocks of the 
propagator have the same structure as in the case of the 
2SC/g2SC phase. 

From Fig. we see that there is indeed a gapless 
mode in the "green-strange-blue-down" sector. This is 
the same that was found in Ref. ^5|- Note also that, 
in agreement with Ref. fl5||. the "red-strange-blue-up" 
quasiparticle has a dispersion relation that is nearly 
quadratic, k rb u ~ \k - k*\ 2 with k* m 400 MeV for a 
given choice of parameters, see Fig. (c). The nearly 
quadratic dispersion relation resembles the situation at 
the transition between the 2SC phase, where n r u = n d , 
and the gapless 2SC phase, where n r u and n d are diffe- 
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rent. This explains the approximate equality n r a « n b u 
mentioned above. 



IV. RESULTS AT NONZERO TEMPERATURE 

In this section, we present the results for the phase 
structure of dense neutral three-flavor quark matter in 
the plane of temperature and m^/jit, as well as in the 
plane of temperature and quark chemical potential. 

Let us start with the discussion of the temperature de- 
pendence of the gap parameters in the two qualitatively 
different cases of small and large values of the strange 
quark mass. As we saw in the previous section, the zero- 
temperature properties of neutral quark matter were very 
different in these two limits. 

The results for the temperature dependence of the gap 
parameters are shown in Figs. and El for two diffe- 
rent values of the strange quark mass that represent the 
two qualitatively different regimes. In the case of a small 
strange quark mass (i.e., the case of rr? s j[i = 80 MeV 
which is shown in Figs. and the zero-temperature 
limit corresponds to the CFL phase. This is seen from the 
fact that the three different gaps shown in every panel 
of Fig. [3 merge as T — > 0. At nonzero temperature, on 
the other hand, the gap parameters are not the same. 
This suggests that, similar to the zero-temperature case 
of Figs. and a phase transition to the gCFL phase 
happens at some nonzero temperature. However, we shall 
show below that there is no phase transition between 
the CFL and gCFL phases at any nonzero temperature. 
Instead, there is an insulator-metal crossover transition 
between the CFL phase and a so-called metallic CFL 
(mCFL) phase. At this point, all quasiparticles are still 
gapped. At some higher temperature, the mCFL phase 
is replaced by the gCFL phase. 

If the temperature is increased even further, there are 
three consecutive phase transitions. These correspond to 
the three phase transitions predicted in Ref. |16| in the 
limit of a small strange quark mass. In order to resolve 
these, we show a close-up of the near-critical region of 
Fig.0in Fig.0 The three transitions that we observe are 
the following: (i) transition from the gCFL phase to the 
so-called uSC phase; (ii) transition from the uSC phase 
to the 2SC phase; (ii) transition from the 2SC phase to 
the normal quark phase. Here, the notation uSC (dSC) 
stands for superconductivity in which all three colors of 
up (down) quark flavor participate in diquark pairing 
[l6j . Our results differ from those of Ref. in that 
the dSC phase is replaced by the uSC phase. The rea- 
son is that, in our case, the first gaps which vanish with 
increasing temperature are 4>\ and tpi , see Fig. while 
in their case A2 (corresponding to our (f>2 and ^2) dis- 
appears first. Although we use the same terms for the 
phases that were introduced in Ref. llj, we distinguish 
between the gapped phases (e.g., CFL and mCFL phase) 
and the gapless phases (e.g., gCFL). Also, in order to 
reflect the physical properties of the mCFL phase, we 



suggest to use the term metallic CFL, instead of modi- 
fied CFL as in Ref. [l^. (Note that, in that work, the 
mCFL phase also encompasses the gCFL phase.) 

In the case of a large strange quark mass (i.e., the 
case of m^/fi = 320 MeV shown in Fig. 01, the zero- 
temperature limit corresponds to the gCFL phase. By 
looking at the corresponding temperature dependence of 
the gap parameters, we see that this case is a natural gen- 
eralization of the previous limit of a small strange quark 
mass. There are also three consecutive phase transitions. 
It is noticeable, however, that the separation between the 
different transitions becomes much wider at large m s . 

We now take a closer look at the transition between 
the CFL, the mCFL, and the gCFL phase. Let us recall 
that, at zero temperature, there was no symmetry con- 
nected with the order parameter, i.e., the number density 
of electrons, that is associated with the CFL — > gCFL 
phase transition. At nonzero temperature, the electron 
density is not strictly zero in the CFL phase as soon as 
m s 7^ 0. Indeed, the arguments of Ref. |3Ij regarding the 
enforced neutrality of the CFL phase do not apply at 
T^O. This leads us to the conclusion that the insulator- 
metal transition between the CFL and the mCFL phase 
is just a smooth crossover at T 7^ 0. Of course, in princi- 
ple, we can never exclude the existence of a first-order 
phase transition. Our numerical analysis, however, re- 
veals a crossover. The transition can only be identified 
by a rapid increase of the electron density in a relatively 
narrow window of temperatures, see Fig.^|(a). The lo- 
cation of the maximum of the corresponding "susceptibil- 
ity" (i.e., driQ/dT) is then associated with the transition 
point. 

The transition between the mCFL and the gCFL phase 
corresponds to the appearance of gapless quasiparticle 
modes in the spectrum. We do not yet know whether this 
transition is associated with any physical susceptibility. 
There is no way of telling from the temperature depen- 
dences in Figs. and whether the corresponding 
CFL and/or 2SC phases are gapless or not. This addi- 
tional piece of information can only be extracted from 
the behavior of the quasiparticle spectra. We also inves- 
tigated them, but we do not show them explicitly. 

Our results for the phase structure of dense neutral 
three-flavor quark matter are summarized in Fig. If II We 
show the phase diagram in the T-m 2 s / \i plane at a fixed 
value of the quark chemical potential, [i = 500 MeV, and 
in the T-fj, plane at a fixed value of the strange quark 
mass, m s — 250 MeV. The three solid lines denote the 
three phase transitions discussed above. The two dashed 
lines mark the appearance of gapless modes in the mCFL 
and 2SC phases. We could term these as the mCFL — > 
gCFL and 2SC — > g2SC crossover transitions. In addition, 
as we mentioned above, there is also an insulator-metal 
type transition between the CFL and mCFL phase. This 
is marked by the dotted lines in Fig. ^2 
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V. CONCLUSIONS 

In this paper, we studied neutral three-flavor quark 
matter at large baryon densities. We obtained a very rich 
phase structure by varying the strange quark mass, the 
quark chemical potential, and the temperature. 

At T = 0, there are two main possibilities for the 
strange quark matter ground state: the CFL and gCFL 
phases. These findings confirm the results of Ref. 0] 
concerning the existence of the gapless CFL phase, the 
estimate of the critical value of the strange quark mass 
m s , and the dependence of the chemical potentials on 
m s . We also confirm that it is the color neutrality condi- 
tion, controlled by the color chemical potential which 
drives the transition from the CFL to the gCFL phase 
[l5| . This is in contrast to gapless 2SC superconductivity 
which results from electrical neutrality [llj. 

Because we use a nine-parameter ansatz for the gap 
function, the results of this paper are more general than 
those of Ref. 15]. For example, we were able to ex- 
plicitly study the effects of the symmetric pairing chan- 
nel, described by the sextet-sextet gap parameters, that 
were neglected in Ref. 0. As one might have expected, 
these latter modify the quasiparticle dispersion relations 
only slightly. This check was important, however, to see 
that the zero-temperature phase transition from the CFL 
phase to the gapless CFL phase, which is not associated 
with any symmetry, is robust against such a deformation 
of the quark system. 

In this paper, we also studied the temperature de- 
pendence of the gap parameters and the quasiparticle 
spectra. In particular, this study revealed that there ex- 
ist several different phases of neutral three-flavor quark 
matter that have been predicted in the framework of 
the Ginzburg-Landau-type effective theory in Ref. 
Our results extend the near-critical behavior discussed 



in Ref. |l6( to all temperatures. Also, we show how this 
behavior evolves with changing the value of the strange 
quark mass. The only real qualitative difference between 
our results and the results of Ref. 0] is that, instead 
of the dSC phase, we find the uSC phase in the phase 
diagram. 

The main result of our paper is the complete phase dia- 
gram of neutral three-flavor quark matter in the T-m 2 s //i 
and T—jj, plane, shown in Fig. 1111 In this figure, all sym- 
metry related phase transitions are denoted by solid lines. 
(The symmetries of all phases appearing in this figure 
were discussed in Ref. |l|j-) I n the mean-field approxima- 
tion used here, all of these transitions are second-order 
phase transitions. After taking into account various types 
of fluctuations, the nature of some of them may change 
[32|. A detailed study of this issue is, however, outside 
the scope of this paper. The dashed lines in Fig.lTTlsepa- 
rate the mCFL and regular 2SC phases from the gapless 
CFL and gapless 2SC phases. These cannot be real phase 
transitions, but are at most smooth crossovers. At T = 0, 
there is an insulator-metal phase transition between the 
CFL and the gCFL phase [Hi- At nonzero temperature, 
there exists a similar insulator-metal type transition be- 
tween the CFL and the mCFL phase, given by the dotted 
lines in Fig. HP 
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FIG. 2: The absolute values of the gap parameters as a function of ml/ 'fx for neutral color-superconducting quark matter at 
T — and /i = 500 MeV. The actual values of the gap parameters shown in panel a) are negative. 
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FIG. 3: The electrical and color chemical potentials as a function of m^/fx of electrical and color neutral color-superconducting 
quark matter at T = and fj, = 500 MeV. 
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FIG. 5: The number densities of each quark color and flavor as a function of /x for electric and color neutral color-superconducting 
quark matter at T = and m 3 — 300 MeV. The densities are given in units of the saturation density of nuclear matter, 
po = 0.15 fm" 3 . 
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